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Abstract 

We apply a small magnetic field to strongly interacting matter with a gravity dual 
description as an electron star. These systems are both metallic and quantum critical 
at low energies. The resulting quantum oscillations are shown to be of the Kosevich- 
Lifshitz form characteristic of Fermi liquid theory. It is seen that only fermions at a 
single radius in the electron star contribute to the oscillations. We proceed to show 
that the Fermi surface area extracted from the quantum oscillations does not obey the 
simplest statement of the Luttinger theorem, that is, it is not universally proportional to 
the total charge density. It follows that our system is a non- Fermi liquid that nonetheless 
exhibits Kosevich-Lifshitz quantum oscillations. We explain how the Luttinger count is 
recovered via a field theoretic description involving a continuum of 'smeared' fermionic 
excitations. 



1 Introduction 



Fermi liquid theory is a robust and weakly interacting low energy description of a finite 
density of interacting fermions [U [2]. In order to achieve non- Fermi liquid metallic be- 
haviour, the low energy theory must contain additional gapless degrees of freedom beyond 
the particle-hole excitations of the Fermi surface. Gapless bosonic modes can arise either 
when a Fermi liquid is tuned to a quantum critical point [31 01 [5] or in critical phases wherein 
kinematic constraints require the introduction of gauge symmetries and the associated pro- 
tected gapless excitations [6j [T] . 

In non- Fermi liquids, the gapless bosonic modes generate relevant perturbations of Fermi 
liquid theory. The system therefore flows to a strongly interacting low energy regime. It 
has recently been appreciated that expansions in a large number ./V of fermion species fail 
to halt the flow to strong coupling [8j El [10]. Cousins of the usual (4 — e)-type expansions 
lead to controlled weakly interacting non- Fermi liquid fixed points [TTJ [p2J H3] . In this letter 
we aim to gain insight from an approach to metallic quantum criticality that is not built 
around weakly interacting notions of single particle Green's functions and self energies. 

The holographic correspondence provides a classical dual description of certain strongly 
interacting large N gauge theories. The most established examples are maximally supersym- 
metric Yang-Mills theories [H]. More generally, one expects gravity duals for a 'landscape' 
of field theories [15] . Known examples typically have SU (N) gauge fields coupled to fermion 
and bosonic matter that can in addition be charged under global symmetries. In these cases 
we can place the theory at a nonzero chemical potential fi for a global U{\) symmetry and 
flow to energy scales E -C fj, |16j . If the low energy theory contains a charged bosonic 
operator with low enough scaling dimension, then a superfluid phase transition will result 
|17 } I18| [15]. Alternatively, if the theory contains low lying charged fermionic operators, the 
'bulk' gravitational description is given by a fermion fluid that drives the theory to an IR 
fixed point characterised by a finite dynamical critical exponent z \19\ 120] PI 




A self gravitating fluid of charged fermions in asymptotically Anti-de Sitter space was 

1 If neither sufficiently low lying bosons nor fermions are present in the operator spectrum, the bulk 
is an extremal black hole. Extremal black holes have an AdS2 x R 2 near horizon geometry and therefore 
correspond to z — oo fixed points. Phenomenologically attractive consequences of such fixed points have 
recently been emphasized [211 122] . The zero temperature entropy density necessarily associated with z = oo 
scaling suggests that this feature is likely an artifact of the large N limit. While the AdS2 regime can be 
pushed to parametrically low energy scales by coupling the bulk Maxwell field to dilaton fields 23, 24, 25 , 
eventually higher derivative terms in the action become important and are expected to stabilise the flow at 
a 2 = oo IR fixed point [26| . 
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termed an 'electron star' in [20] by analogy with the familiar neutron stars of astrophysics 
\27\ I28j . which satisfy very similar equations. At zero temperature the entire charge is 
carried by the fermions. The absence of additional charge-carrying sectors of the theory, 
such as extremal black holes, suggests that electron stars are a promising arena to investigate 
the nature of strongly interacting Fermi surfaces. In this letter we will probe two key 
features of Fermi surfaces: their manifestation as quantum oscillations in a magnetic field 
and the Luttinger theorem relating the area of the Fermi surface to the total charge density 
|29j . One of our main results will be that while the system exhibits the classic Kosevich- 
Lifshitz quantum oscillations characteristic of Fermi liquids, the corresponding Fermi surface 
area does not satisfy the Luttinger theorem. While a direct comparison is premature, it 
is noteworthy that this same apparent dichotomy between Fermi liquid oscillations and 
non-Fermi liquid physics has haunted much recent discussion in the cuprates following the 
seminal experiments of [HU]. The geometry behind the mismatch in our holographic setup 
will become apparent: many of the bulk fermions do not participate in the oscillations 
and therefore disrupt the 'Luttinger count'. In the final section we give a field theoretic 
description of this phenomenon and compare with 'fractionalized Fermi liquids' [311 I32| . A 
sketch of an electron star is given in figure [TJ 




Figure 1: The electron star. The field theory charge density is given by the bulk electric field, 
which is sourced entirely by the fermion fluid. Quantum oscillations are due to fermions 
at a single radius in the star. The quantum critical region emerges due to screening of the 
electric field by the fermion fluid. 
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2 Electron stars distilled 



The electron stars of interest to us will be zero temperature solutions of Einstein-Maxwell 
theory with a negative cosmological constant and coupled to a perfect fluid of unit charge 
relativistic fermions. Throughout this paper we will consider a 3+1 dimensional bulk for 
2+1 dimensional field theories. The equations of motion are 

Rab ~ \g ab R - ^9ab = k 2 (F ac F b c - - A g a bF cd F cd ^ + T ab ^j , (2.1) 

and 

V a F ba = e 2 J b . (2.2) 
As usual the perfect fluid energy momentum tensor and current are 

Tab = {p + p)u a u b + pg ab , J a = au a . (2.3) 

These two quantities must be conserved and the four velocity u normalised so that u 2 = — 1. 
The pressure p, energy density p and charge density a are fields on spacetime. They are 
given in terms of a local equation of state by 

/*Mloc. /"Mloc. 

p = / E g(E)dE , a= / g(E)dE , -p = p - p loc a , (2.4) 

J m J m 

where m is the fermion mass, the local density of states is 

g{E) = fiE^E 2 - m 2 , (2.5) 

with (3 an order one number, e.g. (3 = tt~ 2 for two component spinors, and the local chemical 
potential 

A t 

Woe. = ■ (2.6) 



The electron star has the following metric and Maxwell field 

ds 2 = L 2 (-fdt 2 + gdr 2 + (dx 2 + dy 2 )^j , A = ^hdt . (2.7) 

Here f,g,h are functions of the radial coordinate r. The pressure and energy and charge 

densities are also functions of r. The velocity has nonzero component u = l/(Ly/f). The 

couplings can be scaled out of the equations of motion by writing 

1 1 1 e 4 L 2 k 2 

p=— F ^-p, p=——p, a = —^a, P= — g-/?, m 2 = -^m 2 . (2.8) 
L a k a L z k z eL z K k z e z 

The quantities rh and j5 are the two free dimensionless parameters of the theory. 

The resulting equations of motion for {f,g,h,p,p,a} are written down and solved nu- 
merically in [20]. The key features of the solution are that 
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Outside of the star radius, r < r s , there is no fluid present and the spacetime is that of 
an AdS4 planar Reissner-Nordstrom black hole. Thus the solution is asymptotically 
characterised by a chemical potential ft, charge density Q and energy density E. The 



hats over these dual field theory quantities indicate that rescalings analogous to (2.8) 
have been performed. These appear in the geometry as 

c 2 o r 2 Q 2 c 2 
f= Mr + -%-, 9 = — f , h = {x-rQ. (2.9) 



The constant c keeps track of the normalisation of the time coordinate. Zero temper- 

b 



ature together with scale invariance of the UV theory imply that E = ^fiQ. 



Far inside the star, r 3> r s , the solution acquires a scale-invariant form characterised 
by a dynamical critical exponent z 

1 = — Z ^ 9=—, h =—- ( 2 - 10 ) 

The exponent z can be determined numerically in terms of rh and (3. At a given fixed 
rh the exponent grows unboundedly as /3 — > and monotonically tends to a constant 
as /3 — > oo. For all {to, /3} 
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1 — TO 2 

Furthermore in order for the electron star to exist the fermion must not be too massive. 
Specifically < rh < 1. This also ensures that z > 1. 

For the emergent scaling exponent z to be order one rather than parametrically large 
|19j while simultaneously requiring the classical bulk description to be valid, one needs [20] 

e 2 ~y<l. (2.12) 

Recall that L is the cosmological constant scale, e the Maxwell coupling and k 2 the Newton 



constant. The condition (2.12 ) can be rephrased as the following dual field theoretic relation. 
Two 'central charges' can be defined from the zero temperature current J and energy 
momentum tensor T correlations: (TT) ~ c and {J J) ~ k. As usual holography requires 



c> 1. The additional condition (2.12) becomes the statement that c ~ k . The current J 
is normalised here such that the fermion operator has charge one. 

In between the UV and deep IR geometries is the bulk of the electron star. The star 
is an equilibrated self gravitating system of fermions that carry the entire electric charge. 
In the remainder of this paper we proceed to a dual field theory characterisation of the 
fermionic physics of the intermediate region, geometrically associated with the field theory 
energy scale E ~ fj,. 
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3 Magnetic field probe of electron stars 



A density of fermions leads to robust characteristic features in the magnetic field dependence 
of thermodynamic quantities. Expanding in small momenta about a point feo on the Fermi 
surface of a Fermi liquid leads to the Green's function 

k 2 

G F 1 L (u;,k) = -ioj + v*k ± + ^. (3.1) 

Here uj is the Euclidean frequency measured relative to the chemical potential fi. The 
constants v* and m* are the (locally defined on the Fermi surface) Fermi velocity and 
effective electron mass. The momenta k± and fen are the components of the momentum in 
directions perpendicular and parallel to the Fermi surface at the point feo- The dispersion 



relation (3.1 ) about a Fermi surface leads to de Haas- van Alphen quantum oscillations in the 
magnetic susceptibility. Restricting for simplicity to a circular Fermi surface, the leading 
oscillation in the susceptibility per unit area is given by the Kosevich-Lifshitz formula [33] 

d 2 Q fi 2 m* 2ir 2 m*T/B A F 
Xosc. B 2 S inh(27r2m*r/i?) C ° S • (3 ' 2) 

Here T is the temperature, B the magnetic field and Ap = irkp is the area of the Fermi 
surface. The Fermi momentum fep = m*v*. Thus from a measurement of the susceptibility 



one can determine the two coefficients in the dispersion relation (3.1). 

We will now search for an oscillatory magnetic susceptibility of our electron stars. 
Adding fully fledged magnetic fields and temperature dependence is involved and will not 
be considered here. Instead, we will work in a regime where in our field theory T, y/B -C /x. 
This is typically the experimentally relevant regime. In this limit we can ignore the back- 
reaction of temperature and magnetic field on the spacetime geometry as only the very 
near horizon region will be affected. Thus we keep the radial profiles of the functions 
{/, g, h,p, p, a} satisfying the equations of the previous section and obtained numerically in 
|20| . The crucial consequence of the magnetic field is to collapse the local fermion spectrum 
into Landau levels labelled by t = 0, 1, 2.... While B is taken small, IB must be kept finite. 
The splitting of the spectrum leads to non-analyticities in the magnetic field dependence 
that are then smeared out by the temperature - the competition between these two effects 



leads to the formula (3.2) in a Fermi liquid. 

Thus to leading order at low temperatures and magnetic fields we can use the expression 
for the free energy density from [20] 



r Q = ~Y L ~ 2 



n = --fxQ = I v^a(s)ds. (3.3) 



This formula simply expresses the charge density of the dual field theory as an integral 
over the charge density in the electron star. In the presence of a small magnetic field and 



temperature, the local charge density a in (2.4) must now be written as (we drop from the 
outset the contribution of antiparticles, these do not see the Fermi surface and hence do 
not contribute to quantum oscillations) 

P -Bloc. 



a 



/oo i 
dP Yl i + e (£(<,p)-Woc.)/Tloc ' ^ 3 ' 4 ) 



where E(£,p) = \/p 2 + m 2 + 2£B\ oc . . We have dropped a factor of \ that appears in 
the £ = term. The local temperature and magnetic field are related to the rescaled 
temperature and magnetic field of the dual field theory via 

Aoc(r) = = - 2 r 2 B , T loc .(r) = = 6 ^ . (3.5) 

The first of these relations used the fact that in the absence of magnetic sources the magnetic 
flux F xy is constant. In the next few equations we drop the 'loc' subscript on fi, T and B 
for clarity of expression. 

A finite temperature affects electron stars in two ways. It leads to Fermi-Dirac terms 



in the sum over states as in (3.4) and also results in the appearance of a black hole in 
the interior of the geometry. The black hole is classical in the bulk and thereby captures 
the leading finite temperature physics whereas the effects of Hawking radiation, which is 
responsible for exciting the fermion fluid, are suppressed in the classical limit |34j . Changes 
in the local geometry, however, are insufficient to smooth out the quantum oscillations and 



therefore the expression (3.4) contains the leading effects of temperature for our purposes. 



While performing the sums in (3.4) is a standard computation, we will outline the 



method used in |35j which extracts the oscillating part efficiently and transparently. We 
first perform a Poisson resummation of the sum over Landau levels 



a 

2 ' 

k 



H i M d P 1 T AE(L P )-,)/T ■ ( 3 - 6 ) 



1 + e {E(L P )-p,)/T ■ 

We then proceed to change integration variables from £ to the energy E and also to decom- 



pose the Fermi-Dirac distribution into a sum over Matsubara frequencies (the latter involves 
adding a /U-independent term that will be cancelled by the antiparticle contribution). Thus 

RT r+oo poo Trik(E 2 —p 2 —m 2 )/B 

a = -^-Y / dp / dEE —. ^. (3.7) 



2 - -Vp 2 + m2 E — fi (l -\- i— 27r (n 1 1 



k,n "Vf-vm- Cj - fj, -r -t- 2, 

The integral over E is performed by deforming the integration contour in such a way that 
the argument of the exponential is always real and decreasing at large energies. The only 
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obstruction to doing this is the pole coming from the Matsubara frequencies. This obstruc- 
tion is only present if p 2 + m 2 < /i 2 and furthermore n + \ > for > or ti + ^ <0 for 
k < 0. The rotated contour integral is manifestly non oscillating in B. The only oscillating 



contribution comes from the poles and is given to leading order in — by 



(To 



2vr/3T/iIm e 

k,n>0 



nik(^ 2 -m 2 )/ B -4TT 2 k(n-^)Tn/B 



■v/ [i 2 —m 2 



dpe 



-itikp 2 /B 



(3i 



The expression (3.8) can now be plugged into the formula (3.4) for the free energy of 

^2 



the boundary field theory. So long as the local Fermi surface area proportional to /j 



m 



is not small, i.e. so long as we are not at the boundary of the electron star, then the fact 
that fi 2 /B <C 1 allows us to perform a saddle point evaluation of the integral over radial 
position s in (3.4). Because T <C /J, the imaginary exponent in (3.8) in the more rapidly 



varying one. After taking the saddle point we restrict to the first harmonic k = 1 as higher 
modes are suppressed and then perform the sum over the Matsubara frequencies n. Because 



m ^> B, the p integral in (3.8) can be performed over the whole real line. The final 



result for the oscillatory free energy density is the Kosevich-Lifshitz form 

cfiB 2 



a 



>'* 



3V2vr 



2-n z rh*T/B A F 

^=^= 7i — — cos 

fti i r \\ sinh(27r 2 m*T/5) B 



(3.9) 



where we have expressed everything in terms of rescaled boundary field theory quantities. 

irface a 

ir(h 2 /f — nr 



Comparing with (3.2), the effective Fermi surface area is 

,2\ 



A F 



max ■ 



(3.10) 



with TV being the radius at which the maximum is attained, and the effective fermion mass 

h 



m 



fr 2 



(3.11) 



The maximum required in (3.10) exists and is unique. This is because the numerator, the 
local Fermi surface area, increases monotonically from zero at the boundary of the star 
to a constant in the interior, while the denominator is constant on the star boundary and 
diverges towards the interior r — > oo. The physical reason the denominator appears is that 
the conserved magnetic flux is fixed at the boundary of spacetime and becomes concentrated 
towards the interior according to (3.5). The radius r± is shown in figure [2] below for various 
values of the parameters {m,z}, using the numerical methods of [20]. The fact that low 
energy physics is captured at a bulk radius of order \i rather than r™ 1 <C /J may indicate 
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Figure 2: Radius contributing to quantum oscillations over the radius of the electron star. 
Larger r+/r s > 1 is the interior of the star. From top to bottom, curves have rh = 0.07, 0.36 
and 0.7. We are characterising the solutions by their IR scaling exponent z rather than f3. 



that the usual holographic IR/UV correspondence is complicated in the presence of a Fermi 
surface, analogously to the modified RG flow in field theory [TJ [2]. 

The previous three equations are the first important result in this paper. We have shown 
that upon applying a small magnetic field at low temperatures, our strongly interacting field 



theory at finite density exhibits quantum oscillations with a Kosevich-Lifshitz form (3.9). 
Furthermore, we found that the Fermi surface area and effective quasiparticle mass defined 
through these oscillations are determined by the properties of the electron star at a unique 
radius r*. The fact that only an extremal orbit in the Fermi surface contributes is of course 
familar [36jJ^] The novelty in our setup is that the surface is extremal in real space rather 
than momentum space and that the spatial direction in question is not part of the field 
theory itself, but rather the emergent dimension characteristic of holographic duality. The 
general relativistic physics of the extra dimension manifests itself in e.g. the appearance of 



the metric function / in (3.11 ): The effective mass and Fermi surface area are not related to 
each other as for free electrons with a chemical potential fx, there is a gravitational redshift 
in the effective mass due to the fact that it is defined at a fixed radius in the bulk. 

The dependence of Ap and rh* on the parameters of the model {rh, z} can be computed 
using the numerical techniques of [20] and are shown in figure [3] below. It is important 



2 This is seen clearly at zero temperature: A nonanalyticity in the free energy arises whenever a Landau 
level goes from being occupied at some radii to being above the local chemical potential at all radii. The 



nonanalyticity is therefore localised at the extremal radius (3.101 at which the Landau level is last occupied. 



to plot dimensionless quantities and rescale by factors of c in (2.9) to ensure that theories 
are compared with the same time coordinate. The left hand curve shows the Fermi surface 




Figure 3: The 'Luttinger ratio' of the Fermi surface area to charge density (left) and the 
effective quasiparticle mass over the chemical potential (right) as a function of IR scaling 
exponent z. From left to right, curves have rh = 0.07,0.36 and 0.7. 



area divided by the charge density. This is a dimensionless quantity without any factors 
of the speed of light. Restoring couplings Ap/Q = ^^Ap/Q <C 1. Although somewhat 
dependent on the convention for charge normalisation, this is consistent with the general 
holographic intuition that the field theory is in a 'large N' limit, so that there many species 
of fermion and hence a large charge density with a fixed Fermi surface area. We call this 
quantity the 'Luttinger ratio' as for a Fermi liquid the Luttinger theorem states that it is 
constant independently of the strength of interactions in the microscopic theory [29J. The 
important conclusion of the left plot in figure [3] is that the Luttinger theorem does not 
hold for these systems, despite the fact that the quantum oscillations obey a Kosevich- 
Lifshitz formj^] In the plot we have varied the UV dimension of the charged operator, via 
rh, and a ratio of the current and energy momentum central charges, via j3, see below (2.12) 
above, keeping the charge of the fermionic operator constant. One might worry whether 
the normalisation of the charge operator should be changed simultaneously in a such a way 
that the Luttinger theorem holds. However, we will now explain the physics behind the 
violation of the Luttinger theorem in electron stars and their dual field theories. 



3 Previous holographic considerations of the Luttinger theorem appear in the context of probe fermions 
in an extremal black hole background |37l 138] . Without backreaction, the entire charge is carried by the 
black hole and the multiple probe Fermi surfaces include cases for which the quantum oscillations are not 
Kosevich-Lifshitz [391 135] . It is therefore not clear why it is appropriate to ask for a Luttinger theorem. 
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4 Field theory: Smeared Fermi surface and Landau bands 



From the bulk perspective it is clear what is going on. In the bulk we have a Fermi surface 
in a 3+1 dimensional spacetime. Because the spacetime is not homogeneous, the Fermi 
surface is not uniform in position space. As is well known [36J and as we noted above, 
quantum oscillations are only sensitive to extremal cross sections of the Fermi surface. If a 
3+1 dimensional Fermi surface is not simply a sphere in momentum space, then the area of 
the Fermi surface measured by quantum oscillations does not contain enough information 
to deduce the full volume of the Fermi surface. This is manifestly what is happening in our 
bulk setup: only fermions at an extremal radius are contributing to quantum oscillations. 
Hence the resulting Fermi surface area cannot determine the full volume of the bulk Fermi 
surface, which would need to be integrated over all radii. In our fluid limit the Luttinger 
theorem does hold locally at each radius. Therefore integrating the total Fermi volume 
over all radii would indeed reproduce the total charge density as required by the Luttinger 
theorem. 

From the perspective of the dual strongly interacting field theory the interpretation is 
much more interesting. The dual field theory is 2+1 dimensional and spatially homogeneous. 
Therefore any Fermi surface will be a circle in momentum space and quantum oscillations 
should provide a unique definition of the area of this Fermi surface. The question becomes, 
in the spirit of [29], what low energy degrees of freedom are not captured by the quantum 
oscillations and are thereby responsible for violating the 'Luttinger count'? 

The bulk fermion field we have discussed so far corresponds to a 'single trace' fermionic 
operator in the dual field theory. In our 'large iV' limit, single trace operators have 
correlators that factorise despite the fact that they do not obey linear wave equations. 
Thus if we decompose the operator into energy and momentum modes 

= / ^ + c^e*-*-*-) , (4.1) 

then we can expect a useful description in terms of the creation operators , . Essentially 
the same decomposition is considered in [40, 41]. Whereas those papers work on a spatial 
sphere and use the state-operator correspondence, we are using the existence of a Fock space 
at large N. Unlike for free fields, however, we should not expect that only operators with a 
simple field theoretic 'on shell' relation between uj and k contribute to physical quantities. 
Rather the on shell condition is that of the higher dimensional bulk theory 

g tt (r)u 2 + g xx (r)k 2 + g rr {r)k 2 r = m 2 . (4.2) 
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This shows that for a given field theory momentum k, fermions with a range of u will be 
important. We can subsume the two bulk label^Av and r into a single label, the effective 
mass M of the field theoretic fermion. On shell fermions will then have 

u) k (M) = ^c(M) 2 k 2 + M 2 . (4.3) 

In the limit in which we are working, bulk interactions are such that fermions modify the 



background geometry, but the fermions still satisfy the Gaussian dispersion (4.2). 

We can now imagine integrating out all degrees of freedom except for the charge-carrying 
fermionic field The effective action that captures the fermionic ground state of the theory 
is then given by a continuum of modes labelled by their mass M 

H cS . = J dMA(M) J d 2 k (w fc (M) - //) c{(M)c k (M) , (4.4) 

where in terms of the modes (4.1) of the operator 

c fe (M) = c Wfc(M)>fc . (4.5) 

The functions A{M) and c(M) in the action are determined in the process of integrating out 
the other fields. In the dual bulk description this amounts to solving the Einstein-Maxwell- 
fluid equations as was done in [20] . The functions A(M) and c(M) can then be read off 
from the distribution of the fermion fluid in the bulk radial direction and the background 
spacetime geometry. The effects of the emergent IR scaling are contained in these functions. 
Because the single operator \E r does not lead to a single fermionic excitation but rather a 



continuum thereof, we might refer to the c k (M) excitations of (4.4) as 'smeared' fermions. 



Given the action (4.4), it is clear that the Luttinger count and quantum oscillations are 
reconciled in the same way as in the bulk description. The Fermi surface is labelled by both 
M and the momentum. Only fermions that are at an extremal cross section in this total 
space contribute to quantum oscillations. Said differently, there is a continuum of Landau 
levels with an energy depending on M. These form 'Landau bands'. Oscillations are due 
only to electrons at the extrema of the Landau bands. All of the fermions however contribute 
to the total charge density. Previous methods of violating the Luttinger theorem involve 
topological order, in which fractionalised excitations carry off part of the UV electronic 
charge density [3H [32] . Integrating out these excitations may lead to a structure similar 



to the effective smeared action (4.4). It would be very interesting to develop parallels. 
The holographic setup has the virtue that the smeared spectral weight can be described by 
(dually) free stable excitations. 



4 The presence of two labels is due to the coarse grained nature of the bulk ideal fluid description. 
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The effective action (4.4) for the fermions does not capture all of the interesting low 
energy physics. In particular, the dynamical critical exponent z in the IR can be thought 
of as being due to Landau damping of the massless gauge fields of the UV quantum field 
theory by the ground state density of fermions. We will elucidate this connection elsewhere. 
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